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We show that the matrix element of a local quark-gluon operator in the 
■ photon state, (j(kX)\0\'y(kX)) , can be calculated in lattice QCD. The result 

CN . is generalized to other quantities involving space- like photons, including the 

transition form factor 77* — > 7r° and the virtual-photon-nucleon Compton am- 
£> \ plitude (7*iV|7*./V) which can be used to define the generalized Drell-Hearn- 

Gerasimov and Bjorken sum rules. 
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Ever since the 1960's, it has been well known that the photon is not just a point-like 
particle; rather, it has a complicated hadronic structure. Similarity of the photon-nucleon 
scattering cross section to that of meson-nucleon scattering has led to the so-called vector 
dominance model of photon-hadron interactions ffl]. In recent years, the photon structure 
has been explored extensively in e + e~ and ep collider experiments @. For instance, in the 



certain kinematic region, e e scattering is effectively the real and deeply-virtual photon 



collision, from which the partonic structure of the real photon can be extracted. Indeed, 
the unpolarized quark and gluon distributions in the photon has been phenomenologically 
determined to a reasonable accuracy 0, and future data from the polarized HERA and 
e + e~ colliders can constrain the polarized distributions as well B. 

Given the large amount of experimental data available to describe the hadronic properties 
of the photon, it is appropriate to ask how they can be understood from the fundamental 
theory — quantum chromodynamics (QCD). To the authors' knowledge, the answer is still 
open. Part of the reason is that the photon is not an eigenstate of QCD; the standard 
lattice QCD method used in calculating the matrix elements of the nucleon and pion does 
not apply |4[|. As we shall explain, the QCD matrix elements in the photon state are time- 
dependent correlations which are notoriously difficult to access through the Euclidean space 
||. However, we will show in this paper that they can be extracted from the Euclidean 
correlation functions on a lattice. In particular, the moments of polarized (unpolarized) 
quark and gluon distributions of the photon can be calculated in Monte Carlo simulations 
just like for those of the nucleon. Generalizing our finding, we show that a number of 



1 



processes involving space-like virtual photons, such as transition 77* — > 7r and the virtual- 
photon-nucleon forward Compton scattering j*N — > 7*iV, can be studied in lattice QCD. 
The amplitude for the latter process is the key for generalizing the well-known Drell-Hearn- 
Gerasimov and Bjorken sum rules for the spin-dependent structure function of the nucleon 

0- 

We start with a simple but instructive discussion about what is meant by the hadronic or 
QCD structure of the photon. If O(0) is a quark-gluon operator at the spacetime point 0, its 
matrix element in the photon state can be defined from the standard Lehmann-Symanzik- 
Zimmermann (LSZ) reduction formula, 

(7(fcA')|O(0)|7(fcA)) = - lim / d 4 xd 4 ye- lk,( - y - x) k' 2 (0\Te*(\') ■ A(x) 6(0) e(A) • A(y)\0)k' 2 , 

(1) 

where fc M is the on-shell (k 2 = 0) photon momentum, A and A' are the photon helicities, 
and k'^ on the right-hand side is taken to be on-shell after cancelling the photon poles. 
The normalization for the photon states are taken to be covariant, i.e., (7(/c'A')|7(A;A)) = 
2/c (27r) 3 5(fc' — k)5\\/. The renormalization constant Z 3 = 1 + 0(a em ) has been omitted for 
simplicity. Evaluating the electromagnetic part of the Green's function in the lowest-order 
perturbation theory (this is allowed because of the small electromagnetic coupling), we find 

( 7 (fcA')|d(0)| 7 (*A)) = -e 2 J dW 2 /e-^)e^(A')(0|TJ M (o;)6(0)J,( 2 /)|0) e l '(A) , (2) 

where e is the renormalized electric charge unit and is the electromagnetic current op- 
erator of the quarks. According to the above equation, the matrix element of interest is a 
Minkowski-space correlation function in QCD. This is in contrast to an analogous nucleon 
matrix element, (p\0\p), which has no explicit reference to the Minkowski time. 

One can remove the explicit time dependence by deriving a Lehmann representation of 
Eq. (0): separating different time-orderings, inserting complete sets of hadronic states, and 
integrating out the spacetime coordinates x and y. The result can be shown as the time- 
independent perturbation diagrams in Fig. 1. Diagrams a) to f) represent contributions from 
six possible time-orderings. Diagrams g) to i) are the vacuum contributions which must be 
subtracted. As an example, the diagram a) corresponds to the following time-independent 
matrix element, 

Mh\>Mn(ri\MlrH\ r-gy (0|£ * (A/) • Ak)\n')(n'\d(0)\n)(nHX) ■ J(k)\Q) 

(7(fcA J|O(0j|7(fcAj) Fig . la =e ^ — — — , (3) 

where \n) are the eigenstates of Hqqd with momentum k, and to = \k\ is the on-shell photon 
energy. The contributions from the other diagrams can be expressed similarly using the 
rules from time-independent perturbative theory. 
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FIG. 1. Time-ordered diagrams for the quark-gluon matrix elements of the photon state. The 
circles with crosses represent local quark and gluon operators, and the open ones the electromag- 
netic currents. Multiple parallel lines indicate the intermediate hadronic eigenstates. 

Though hardly new, the Lehmann representation helps us in understanding the physical 
origin of the vector dominance. With certain probability, the photon is a superposition of 
a set of hadron states with the photon quantum numbers. The lowest-lying such states are 
7T7T scattering states in the p-wave. A strong resonance in this channel is the p meson. The 
vector dominance then refers to the dominance of the p resonance in the above hadron-state 
summations, although other hadron states can and do contribute. It is interesting to note 
from Fig. 1 that the hadronic part of the photon wave function does not consist of just the 
pure hadron states; it contains also the components that are the mixtures of one or more 
resolved photons and hadrons. 

The above discussion makes it clear that the hadronic structure of the photon is qualita- 
tively different from that of ordinary hadrons, at least in the context of QCD. The photon 
matrix elements contain explicit sums over the complete set of QCD eigenstates as well 
as the contributions from different time-orderings. Because the time-dependent correla- 
tion functions are generically difficult to study in the Euclidean space, there are no general 
rules to calculate them in lattice QCD. One might be tempted to use the electromagnetic 
current J M as a photon interpolating field to extract the matrix element of interest from 
(Q\Jn(y)O(Q)J„(x)\0) at large positive —£4,2/4. The actual outcome is a matrix element 
between the threshold two-7r states, or the p if it were stable. 

The key observation for a lattice calculation of the photon matrix element is that the 
photon is an eigenstate of the combined eletromagnetic and strong interaction hamiltonian 
||. When both interactions are present, we can use the electromagnetic potential as 
the interpolating field for photon. Thus we consider the following matrix elements in the 
Euclidean space, 
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(0|^(y)O(0)A,(x)|0> , (4) 

where x and y are Euclidean coordinates. After Fourier-transforming the spatial coordinates 
and letting y 4 and — x 4 be large, we have 

M^(y 4 ,x 4 ,k) = J d 3 xe~ rk - s J d 3 ye l1 * (0\A^y, y 4 )6(0)A u (x, x 4 )\0) 

-> E(0|^(^)l7(fcA'))e-^(7(fcA')|6(0)| 7 (fcA))e^(7(^A)|^(fc)|0) . (5) 

XX' 

The summations over A and A' can be selected by multiplying M M „ with appropriate po- 
larization vectors. In the second line, we have neglected the three-photon states with the 
same momentum which are suppressed by powers of a em , and the hadronic states which are 
suppressed by an energy gap (2m n at k ~ 0). 

Before putting M Mi ,(y 4 , x 4 , k) on a lattice, we integrate out the photon field in perturba- 
tion theory, 

M^(y 4 , x 4 , k) = e 2 f d^wdhe^^ D w {k, y 4 - w 4 ) 

x(0\T E J p {w,w 4 )d{0)J a {z,z 4 )\0)D au {k,z 4 -x 4 ) , (6) 

where D^ v is the photon propagator in the Euclidean space and T# is the Euclidean time- 
ordering. Using the Fourier integral {k > 0), 

we can write 

M fll/ (y 4 , x 4 , k) ~ J ^ 4 ^ 4 e-" l ^- X4| e-" |y4 - u ' 4| (0|T i; J / ,(-fc,«; 4 )()(0)J i ,(fc^4)|0) . (8) 
Taking the limit y 4 , —x 4 — > 00 and comparing the result with Eq. (|), we get, 
( 7 (A;A , )|O(0)|7(A;A)) = e 2 / dw^^^^e^A'XOlTB J^-jfc, 104)0(0) JAk, ^ 4 )|0)e.-(A) , (9) 



where we have restored the photon polarization in physical gauge (e° = and e-k = 0), i and 
j sum over spatial indices, and J M = J2f e f4>fl[i4>f is the electromagnetic current summing 
over all quark flavors /. The above equation is one of the main results of this paper. 

Several comments on Eq. (H) are as follows: First, the Euclidean correlation function 
(0|TeJ m (&;, w 4 )O(0)J u (k, z 4 )\0) can be calculated in the standard Monte Carlo simulations. 
Hence, the photon matrix elements are accessible in lattice QCD just like the nucleon ones. 
Second, the above expression is nothing but a straightforward analytical continuation of 
the expression in Eq. (0). However, the process arriving at Eq. @ makes it clear that 
the analytical continuation fails if there is no energy gap between the photon and hadronic 
states. Finally, (0|TJ M (fc, w 4 )O(0) J v { k, z 4 )\0) is what one would naively consider when using 

as a photon interpolating field. On the other hand, the correct result involves integrations 
over w 4 and z 4 , which include not only the lowest-energy hadron state but all others with the 
photon quantum numbers. Finally, it is straightforward to show that the above Euclidean 
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correlation function reproduces the Lehmann representation of the photon matrix element. 
In the process, one finds that if the energy of the photon exceeds that of the lowest hadronic 
state (the gap vanishes), Eq. (^) is divergent. Thus, we emphasize that the photon matrix 
element is calculable on a lattice only when there is an energy gap between the photon and 
hadron states. 

In the remainder of the paper, we explore under what conditions the above result can 
be generalized. First of all, the initial and final photons can have different 3-momenta, and 
then the forward matrix element is generalized to form factors. Next, we introduce the 
hadronic structure of an off-shell photon for which cu 7^ \k\. The parton distributions in an 
off-shell photon can be extracted from experimental data just as for an on-shell photon [0. 
In QCD, the off-shell photon matrix elements are defined in terms of the same Minkowski 
correlation in Eq. (0). The corresponding Euclidean expression in Eq. (0) is valid only when 

< lo < \J Am 2 + k 2 , where the upper limit refers to the hadron production threshold for 
a time-like photon. Finally, one can define and calculate on a lattice the hadronic matrix 
elements of the weak-interaction gauge bosons, off- or on-shell, by replacing in Eq. (||) the 
electromagnetic currents with the weak-interaction currents. 

To make further progress, we explore under what conditions a Minkowski correlation can 
be obtained from a numerical calculation of the corresponding Euclidean correlation. Eu- 
clidean Green's functions are generally real as they can be expressed as real integrals in the 
functional formulation and amenable to Monte Carlo simulations, whereas Minkowski corre- 
lations can have imaginary parts corresponding to physical intermediate states propagating 
infinite distances. Therefore, it is reasonable to expect that Minkowski correlations can be 
analytically continued and calculated in Euclidean space only when they are in kinematic 
domains where no physical intermediate states are allowed. In the example of the matrix 

elements for an off-shell photon, the condition cu < \J Am 2 + k 2 ensures that the photon can 
never turn into a real hadron state. 

Thus the transition form factor (amplitude) corresponding to 77* — > 7r° is calculable on 
a lattice, where the first photon is real and the second is space-like. In Minkowski space the 
amplitude can be expressed as 

A^(P,k) = J dSe lky (U\TJ^y)JM\<P)) • (10) 

Denoting the interpolating field for the pion as 0, we construct the following Green's function 
in the Euclidean space, 

G^(P,x 4 ,k) = J dy^ u (0\T E J^ t y 4 )J v {0)<f>{P,x A )\0) . (11) 

— * 

Taking x 4 — > — 00 and keeping only the pion state of momentum P, we have 

G^(P,x 4 ,k) = J ^ 4 e^(O|T £ J M (fc, 2/4 )J i ,(O)| 7 r(P))e^(7r(P)|0(P)|O) , (12) 

where Ep = \J P 2 + m 2 is the on-shell pion energy. By working out its Lehmann representa- 
tion, the Euclidean correlation function reproduces the corresponding Minkowski amplitude. 

The final example involving the photons is the forward virtual-photon-nucleon Compton 
scattering amplitude, 
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T^{P,q) = i J e il}X d 4 x(P\TJ^x)J u (0)\P) , (13) 

where |P) M is the nucleon state with momentum P. The space-like virtual photon momentum 
has two distinct kinematic regions, (P + q) 2 > M% and (P + q) 2 < Mjy-, where is 
the nucleon mass. In the former region, the amplitude is complex and its imaginary part 
corresponds to the inclusive production cross section in eP scattering. In the other region, 
the Compton amplitude is real and therefore we expect it is calculable in lattice QCD. Define 
the four-point Euclidean correlation, 

G^{P,y 4 ,z 4 ,$) = J dx^iOlTEviP^M^x^J^fjiP^z^lQ) , (14) 

where rj is an interpolating field for the nucleon. Take the limit y^, —z^ — > 00, we have the 
nucleon state dominance, 

G^(P,y 4 ,^ 4 ,g) = J d£ 4 e^(0|?7(P)|P>e-^ . 

(15) 

Integrating over x 4 , we find that the remaining correlation function reproduces the physical 
Compton amplitude in the region (P + q) 2 < Mjj. 

A particularly interesting application is when oj — 0. The spin-dependent Compton 
amplitude as a function of Q 2 = —q 2 defines an integral over v for the G\(u, Q 2 ) struc- 
ture function of the nucleon. In the limit Q 2 —>■ 0, the amplitude can be obtained by the 
low-energy theorem and chiral perturbation theory and the integral defines the Drell-Hearn- 
Gerasimov sum rule and its generalization 0. In the limit Q 2 — ► 00, the amplitude can be 
studied in perturbative QCD and the integral defines the Bjorken sum rule and its general- 
ization. Our result here allows a lattice QCD calculation of the amplitude at intermediate 
Q 2 where there is no firm theoretical prediction available. 

In conclusion, we have shown that the QCD structure of the photon can be studied on a 
lattice by calculating Euclidean correlation functions involving local quark-gluon operators 
and electromagnetic currents. Also calculable are a number of other interesting examples 
involving photons, such as 77* — > n° and 7*iV — > j*N. These results shall lead to more 
exciting interplay between lattice QCD and strong-interaction phenomenology. 
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